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FOREWORD 


The determination of opacities and equations of state for astrophysical and 
laboratory plasmas depends upon accurate partition functions. Theoretical gas 
characteristics of plasmas of seven different chemical compositions have been calculated 
and are presented in tabular form. The seven cases are six pure gases (H, He, N, O, Ne, 
and Ar) and the solar chemical composition. All data are tabulated for 10 electron 
pressures from 10-2 dyne/cm? to 107 dyne/cm? and for 50 temperatures from 3000 K 
to 150000 K. The tables include electron density, particle density, the number of free 
electrons per atom, the number of negative hydrogen ions per hydrogen atom, mass 
density, gas pressure, the equation of state parameter (P,/eT), the degeneracy parameter, 
the ionization energy per gram, the ionization plus excitation energy per gram, the total 


energy per gram, the partition functions, the ionization fractions, and the depression of 
the ionization potentials. 
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ANALYSIS 
The partition function is given by 
. kT 
BAT.P= DW Ts Pde ine a) 
=I HY 


where the subscripti refers to the chemical element, the subscript/ refers to the jth ionization state of 
element i, and Six is the statistical weight of the kth energy level, whose excitation potential is Xijkee 


The quantity Wi fT, P,) is the probability that the kth energy state of the ion exists (Refs. 1,2) and is 
given by 


, | | 
log Wi,(T,P,)= 14.69 - Flog N,(T, P,) + 4 log Zp, ~ 4 log ny <0 me) 


where N (7, P,) is the number of quasi-static Stark perturbers per em?, Lig is the effective nuclear 
charge seen by the valence electron, and 1, is the principle quantum number of the Ath energy state. 
The sum in Eq. (1) is terminated at k= m, where m is the first value of n, such that Wir < 1073 and 
B, is accurate to better than 1%. In general, where Way, <1.0,n, =k 


The depression of the ionization potential is given by 
-1 y 


k=m mM. 
oe : : 
AE >) )O- WadM ya [Ea- Xp) * Ey- X01) » AC "Me 8 
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where Ey is the ionization potential of the ion and p is defined as the first value of m, such that 
W.,< 1.0. 
ijk 


The relative number of atoms in successive ionization states is given by the Saha equation 


X,= a Fy ex at cag ae 0 
7 Nig-ry : 
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iG-1) 
where 77 is “the degeneracy parameter and E, is a zero point energy shift. The ionization fraction, n,,, 
is then given by N, g/2Nj- At a given temperature and electron pressure, 7 is given by 


axon" 
a’ Omen ary 2 Fs/9(n) (5) 


where F’3,.(n) is the Fermi-Dirac integral of order 3/2. Then the electron density is given by 


(2m, kT)??? 
N,= ss ae aes F y(n) (6) 
where F’, 2) is the Fermi-Dirac integral of order 1/2. The zero point energy shift, E,, is defined such 
that an electron with zero kinetic energy has zero total energy and is given by (Ref. 3) 


rn 3 
SF Ny Ap yKT(S . 0.731424) 
E, =——= == (7) 
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where Ay = Ze?/ 2A, kT and Ay is the Stewart-Pyatt ion radius (Ref. 4). The Stewart-Pyatt ion radius is 
essentially an interpolating radius between the ion-sphere and Debye-Hiickel approximations. 

The quantities 8 iy Xing Hy and Lik are based upon atomic data from Ref. 5 and Mayer’s opacity 
work (Ref. 3). For energy levels edtside the available atomic data tables, the hydrogenic approxima- 
tion is used. The excitation energy is given by 


eee 
FE emma 8 
where Rhe = 13.5977 eV = 1 Ry and Z;, is evaluated with the use of screening constants from the 

Mayer tables (Ref. 3). For large n’s Mayer’s tables are exceeded, but in this case Zing = =j+1,ie., the 


charge associated with the next higher ion. 


ancnainein antiga 


The number of quasi-static Stark perturbers, N(T, P,), is given by the sum of the number of ions 
plus the number of electrons 


NAT, P,) = dN ir Py? wih P) 


2 
where N, iT, P 2) | is evaluated by direct summation over the ionization fractions. 
In order to enhance convergence, when the quantity n,, <107", we approximate 
72 
AE, = ars (10) 


where S$ = 5.594 X 103/3/4(N, +.N,)''®, i.e, the value of n, such that W,,(7,P,) = 1/2 


The general references for this material are Ref. 6 and a forthcoming article by Fischel and Sparks 


RESULTS 


Since 7 is uniquely determined by the independent variables [Eq. (5)] and the electron density is 
uniquely determined by 7 and T [Eq.(6)], there is only one table of these functions which applies to 
all seven cases. The table of the number of negative hydrogen ions per hydrogen atom is the same for 
the pure hydrogen case and the stellar atmosphere case; therefore only one such table is provided. In 
order to make all the tables as legible as possible, most of the tables have been truncated and printed 
in logarithmic form. The partition functions are printed for those temperatures for which the parti- 
tion function at some electron pressure exceeds the low temperature value, which may differ from the 
ground state statistical weight due to low-lying energy levels. A table of ground state statistical weights 
used in these calculations (Table 1) is provided. The calculations were carried out on an IBM 360 
whose floating point arithmetic is restricted to the number range 10°” to 10%. Values of ny < 1078 
are automatically set to zero and the decadic logarithm of such cases is arbitrarily set to -80. Only 

the temperature cases at which log,, Ny > -80 are printed. Also, only those cases of AE,, for which 

1og19 %z41) 7 ~80 are printed. All ionization states have been calculated. However, due to the above 

printing conventions, not all cases have been printed. Some ions have no tables printed, even when 


higher ions have tables, because of closed or near closed electron shells. An entry of asterisksin any 
table indicates that the case could not be computed. 


*Fischel, D., and Sparks, W. M.: Partition Functions in Ionizing Plasmas. To be published 2 Astrophys. J. 


Table 1—Ground state statistical weights. 
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Appendix I 
LIST OF SYMBOLS 


b’ gas parameter with degeneracy, P, /pT 

B; partition function of element i and ion j 

E, average ionization energy (erg/g) 

creex average ionization plus excitation energy (erg/g) 

Eo zero point energy shift for free electrons (eV) ; 

Ey average total energy (thermal plus ion anon plus excitation) (erg/g) 
Ey ionization potential of ion7 to ionj + 1 for elementi 


AE depression of ionization potential E if due to the ion-ion and ion-electron interactions 


F on), F3/.(n) Fermi-Dirac integrals of order 1/2 and 3/2 respectively 


Six statistical weight of the kth energy level of the ion i,j 
h Planck’s constant, 6.252 X 1072? erg-sec 

i subscript referring to the ith element 

] subscript referring to the jth ion 

k subscript referring to the kth energy level 

k Boltzman’s constant, 1.38046 X 107!6 erg/K 

m value of nm, such that Wiy,.< 16 

m, mass of the electron, 9.1084 X 10728 g 

ny number fraction of element i in ion statej 

Ny principal quantum number of the kth energy level 
N, number density of atomic nuclei (em7?) 

N, number density of electrons (cm™3) 


number density of atoms of element i in the jth ionization state (em™?) 


Xijk 


number density of quasi-static Stark perturbers (em73) 
number density of ions (em73) 

number of negative hydrogen ions per hydrogen atom 
electron pressure (dyne/em?) 

gas pressure (dyne/cm?) 

the Rydberg unit, 13.5977 eV 

temperature (K) 

probability that the Ath energy state of the ion exists 
relative number of atoms of elementi in ionization state j 
number of free electrons per atom N,/N, 


effective nuclear charge seen by the valence electron when it is in energy state k of 
the ion i,j. 


degeneracy parameter 
the number, 3.14159. ..| 
mass density (g/cm?) 


the excitation energy of the kth energy state of the ion ij 
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SOLAR ATMOSPHERE 


Chemical .omposition 


0.858879 X 10° 
0.139310 X 10° 
0.4507 X 1073 
0.8202 X 10-4 
0.7833 X 1073 
0.8589 X 1076 
0.4305 X 1073 
0.1713 X 1075 
0.2158 X 1074 
0.1361 X 10°> 
0.2716 X 10-4 
0.6515 X 10° 
0.4305 X 1077 
0.1228 X 10° 
0.3191 X 107> 
0.9413 X 1077 
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